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I. INTRODUCTION 

There is a growing interest in the hadronic content of the photon |I]] as both a test of 
QCD and as a background to precision measurements of electroweak parameters. In fact, 
some photoproduction processes are dominated by the hadronic interactions of photons in 
certain regions of phase space [0,^. Although there have been several theoretical calculations 
of the hadronic content of the photon @-0, at present the experimental data P] are not 
adequate to significantly constrain theory. Recently a number of papers have appeared 
which present methods of extracting the gluonic component of the photon p|j9|JT0|l. In 
particular, the potential of 67 colliders for the study of the hadronic content of the photon 
has been studied by several authors 0,0. (Such colliders, where high energy e"'"e~ colliders 
are converted into ej colliders, are receiving considerable interest for the phenomenology 
that may be studied in ej collisions.) However, these analysis have concentrated on dijet 
production in 67 collisions where it is difficult to isolate the effect of the gluon and quark 
distributions. A good determination of the quark content of the photon is therefore needed 
to test the various models of the photon structure functions and is essential for the extraction 
of the gluon content of the photon. In this paper we present a novel method of measuring the 
heavy quark content of the photon using hard scattered 6-quarks at e^e~ and 67 coUiders. 
We find that by using b tagging to probe quark distributions inside the photon it should be 
possible to differentiate between the various existing sets of photon distribution functions 
using a 500 GeV e~^e~ collider operating in ej mode. 

Recently, there has been some related work on heavy quark production. In the context of 



leptopro duct ion on a proton, in Ref. |TT| Olness and Riemersma pointed out that there are 
advantages to using both a fixed flavor scheme (in which, e.g., the 6-quark is not considered 
to be a constituent) [T^ and a variable flavor scheme (in which, e.g., the 6-quark may 
or may not be considered a constituent depending on energy scale) |TB[. Although these 
works all deal with leptoproduction on a proton, many of the ideas will apply, with suitable 
modiflcations, to leptoproduction on a photon. In Ref. |T^ , heavy quark production in two 



photon processes at an e~^e~ collider is discussed; although some mention is made of the 
heavy quark content of the photon, the bulk of the analysis deals with the production of 
heavy quark pairs via 77, jg and gg fusion. Finally, Laenen and Riemersma |]l5l present 
a calculation of heavy quark production in ej collisions using a fixed flavor scheme; as we 
consider the heavy quark as being a constituent of the photon, our approach should be 



considered to be complementary to that of Ref. |15 



II. CALCULATION 

We are interested in the measurement of /b/^(x, Q^) which contributes to the process 
e + 7— >e + 6+ jet via the subprocess e + b —>■ e + b which can be written 



cr(e7 -^ efejet) = dx fb/'yix,Q'^) a{eb —> eb). 



(1) 



This process is shown in Fig. 1. It is straightforward to calculate the amplitude for the 
subprocess e + b —>■ e + b, and the squared and summed/averaged matrix element is given 
by: 
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In the above equation, sw = sin 9w and cw = cos 9w with 9w being the weak angle; gyj^ 
are the vector and axial vector couplings of the electron and 6-quark to the Z^, respectively, 
chosen such that g"^ = 1/4 (this fixes the remainder of the couplings unambiguously). Also, 
(5e,b are the electric charges of the electron and 6-quark, respectively. 

In some sense the process 67 — > eb jet can be regarded as an approximation to the direct 
(pointlike) process 67 -^ ebb in analogy to the effective W approximation or the Weiszacker- 
Williams approximation. With appropriate kinematic cuts to ensure that we are describing 
the same kinematic region, the pt distributions of the outgoing b in the two approaches 
closely resemble each other, as demonstrated in Fig. 2. In addition, contributions from the 
process eg -^ ebb, with the gluon coming from the hadronic structure of the photon, are 
similarly contained in our resolved photon contribution via the Altarelli-Parisi evolution of 
the distribution functions. When we consider backgrounds we must take care so as not to 
double count these contributions. 

To determine whether the process is viable as a probe of the quark content of the photon 
we must include detector acceptances and consider the possible backgrounds. Fortunately 
the signature for the process is quite distinct, a single b balanced against the beam electron 
and possibly hadronic jet remnants of the photon. For b detection we use a typical LEP 
detector acceptance for 6's of | cos 6*1 < 0.85 and assume a b detection efficiency of 50%. The 
jet remnants will generally go down the beam so when considering backgrounds we can veto 
events with jets detected above some minimum angle. One sees that the crux of the analysis 
is to determine whether the signal is overwhelmed by backgrounds. 

To illustrate the background contributions we take as an example 6-quark production at 
a 500 GeV ej collider where we have folded in the photon spectrum for the backscattered 
laser and use (to be specific) the Duke and Owens photon structure functions 0. The 
general properties of the other cases we will consider are the same. Possible backgrounds 
are shown in Fig. 3. They can be divided into direct processes (Fig. 3a), 



e+ + e" ^ 6 + 6 (3) 

e+ + e" -> 7 + 6 + 6 (4) 

e + -f ^ e + b + b (5) 

e + 7^z/ + 5 + c (6) 



singly resolved (Fig. 3b) 



e + [g]^ ^e + b + b (7) 

e + [^]^ -^u + b + c (8) 



and doubly resolved (Fig. 3c) 



[gl, + [gh ^b + b (9) 

[q], + [Qh ^b + b (10) 

where we used the notation [p] to represent parton p's content in the photon. We include 
detector acceptance cuts of | cos6'efe| < 0.85 for the b, where 9eb is the angle between b and 
the beam direction, the observed electron (or positron) is at least 10° from the beam and 
all particles not to be detected must be within 10° of the beam except for neutrinos. 

We start with the two processes already mentioned above, 67 — ^ ebb (eqn. 5) and 
e[g]^ -^ ebb (eqn. 7). The pointlike process 67 —>■ ebb contains a coUinear divergence arising 
from the internal 6-quark line which is also included in the process 67 —>■ eb jet. This 



divergence must therefore be first subtracted as described in Ref. |]T3|,|T6| before including 
it as a background. Likewise, the process e[g].y -^ ebb, i.e., with the gluon coming from 
the hadronic structure of the photon, is similarly, to some approximation, contained in our 
resolved photon contribution via the Altarelli-Parisi evolution of the distribution functions. 
A similar collinear subtraction must be made if e[g]^ —* ebb is to be considered a correction 
to our result. We find that, after making the required collinear subtraction, both processes 
67 -^ ebb (direct) and eg — > ebb (once resolved) contribute negligibly in the kinematic region 
under study. 



In the remaining background discussions we are concerned with decay chains that might 
be mistaken for the signaL In particular, in the processes listed in eqns. (3)-(10) heavy 
quarks are pair produced which could lead to the situation where the b (or c) decays semilep- 
tonically and only an electron is seen by the detector. Clearly, this is expected to be rather 
unlikely given that all decay products, both leptonic and hadronic, are expected to be 
boosted in the same direction given the relatively large energy of the produced quarks. The 
likelihood of these events is further reduced given the small branching fractions for these 
decay chains. 

Keeping this in mind we start by describing the possible direct process backgrounds. 
The first two backgrounds are only relevant to the e+e~ cases involving Weiszacker- Williams 
photons. The background e^ + e~ -^ b + b is negligible compared to e"^ + e~ ^6 + 6 + 7 
even though it is technically of lower order in coupling. This is because the b and b come 
out back-to-back with high p^ and are therefore very distinct because the decay products 
of the b are boosted along the b direction making it very unlikely that the electron is seen 
in the detector while the hadronic remnants go down the beampipe. The second direct 
background, e"*" + e^ — i> 6 + 6 + 7, will also turn out to be unimportant in the relevant 
kinematic region once the b is allowed to decay and cuts are imposed so that only the e~ 
and b are observed. The third direct background, e + 7^e + 6 + 6, can contribute in two 
ways. In the first, the beam electron is observed in the detector and the decay products 
of the b are not, while in the second the beam electron goes down the beampipe but the 
electron from b decay is seen in the detector. As discussed above, the situation in which 
the beam electron is hard scattered into the detector is approximated by our signal and 
corrections to the signal, due to this process after subtracting the collinear divergence, are 
found to be negligible in this kinematic region (this process was calculated numerically using 
helicity amplitude techniques). The other possibility, where the beam electron goes down 
the beampipe and is not seen in the detector, can be well described using the Weiszacker- 
Williams approximation and the subprocess 7 + 7^6 + 6 (the subprocess cross section can 
be extracted from, e.g., Ref. [jl^)- The final direct process which could contaminate the 



signal is the charged current process, e + 'j^u + b + c Ref. |TB|; here the electron seen in 
the detector comes from the decay of the c. 

We next consider the singly resolved backgrounds. As discussed above, if the beam 
electron is hard scattered into the detector, the first of the singly resolved backgrounds, 
e + [g]'y ^ e + b + b, (again, we use helicity amplitude techniques to calculate this process 
numerically), is approximated by our signal; corrections to our result in the appropriate 
kinematic region, after the necessary collinear subtraction has been made, are found to be 
negligible. It is also possible that the beam electron goes down the beam pipe, in which case 
this background can be well described in the Weiszacker- Williams approximation with the 



subprocess 7 + (7 ^ 6 + 6 (this cross section can also be extracted from Ref. |T^ by carefully 
modifying couplings and color factors). Also, the charged current process e+ [g]^ — » u + b + c 
can also contribute to the background if the c decays leptonically and the resulting electron 



is observed in the detector; with some care, this subprocess can be extracted from Ref. |]T8 

Finally, we consider the doubly resolved backgrounds where the photon from the beam 
electron contributes a parton from it's hadronic structure. Since the parton model of photon 
structure is defined for real photons only, it is appropriate to use the Weiszacker- Williams 
approximation, with the subprocesses [g]^ + [g]^ —^b + b and [q]^ + [q]^ — > 6 + 6, to calculate 
these backgrounds (see Ref. [|l9l for the cross sections). 

The largest of the backgrounds are displayed in Fig. 4. The p^, (of the b quark) distri- 
butions of the two charged current processes have, as expected, rather long tails due to the 
large mass of the exchanged particle. These p^ distributions are, however, several orders of 
magnitude smaller than our signal, and do not appear on our plot for the scale chosen. The 
remaining subprocesses fall off rapidly with p^ of the b quark; although they may produce 
many events at small p^, they are safely negligible for p^ > 40 GeV. We did not include 
charge identification of the leptons since we include both b and b production in our distribu- 
tions and estimates of event numbers as they both make equal contributions. This approach 
also eliminates the complexities of B^ — B^ mixing in the analysis. 



III. RESULTS 

Having convinced the reader that the signal we are studying is distinct and clean of 
backgrounds we proceed to examine the sensitivity of various kinematic distributions to 
different structure function parameterizations for a number of collider possibilities. We 
consider a group of the existing sets of photonic parton distributions that appear in the 
literature, namely the set of Duke and Owen (DO) IQ, the set of Drees and Grassie (DG) 0, 
the set 1 of Abramowicz, Charchula and Levy (LAC) and the leading order set of Gliick, 
Reyaand Vogt (GRV) 0. 

We first considered the LEP e^e~ collider at CERN upgraded in luminosity using the 
Weizsacker- Williams photon distributions with the 67 cross sections to obtain numerical 
results. The cross section for the signal is expected to be between 0.108 pb (GRV) to 
0.220 pb (LAC). An order of magnitude improvement in the LEP luminosity, resulting in 
1 fm~^/year, would yield 108-220 6's per year or of order 50-100 reconstructed 6's once 
efficiencies are included. This small number of events is unlikely to offer an improvement 
over existing estimates of the b content of the photon. At the LEP200 e+e^ collider the cross 
section is expected to be 0.047 pb (GRV) to 0.092 pb (LAC) which would yield roughly 40 
6's or 20 reconstructed 6's for the total integrated luminosity of 500 pb^^. Clearly a higher 
energy and higher luminosity collider is needed. 

We therefore turn to the proposed higher energy NLC e"'"e~ collider in both e~^e~ mode 
and 67 mode with y^=500 GeV and very high luminosity yielding of order 50 fb^^/year. 
In the e~^e~ mode we fold in the Weizsacker- Williams effective photon distribution and the 
67 mode we fold in the backscattered laser photon spectrum. 

First, consider the e^e~ configuration. The cross section is expected to be between 14 fb 
(GRV) and 24 fb (LAC). Given 50 fb^Vyi' of luminosity, we expect 700 to 1200 events per 
year, and after the 50% b reconstruction efficiency, of order 500 reconstructed 6's per year. 
This will allow a measurement of the 6-quark distribution in the photon at some level, but 
even so, more events are desirable. The Pj, distribution of the 6-quark is shown in Fig. 5a 
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for various sets of photon distribution functions. 

Lastly, consider a y/s = 500 GeV linear e+e^ collider configured as an 67 collider through 
the use of a backscattered laser. In this situation, we expect a cross section of between 68 fb 
(GRV) and 138 fb (LAC), with all backgrounds smaller by at least 2 orders of magnitude 
after the implementation of the cuts described above. This gives a total of 3400 to 6900 
events per year, and after the b reconstruction efficiency, there will still be 2000-3000 events 
per year to analyze. The p^ distribution for the signal is shown in Fig. 5b for various choices 
of photon distributions. 

Having established that only the two cases involving an NLC can provide a useful mea- 
surement of the 6-quark content of the photon {i.e., large enough signal combined with 
manageable backgrounds), we consider useful experimental measurements. 

Before proceeding we point out a subtlety that should be mentioned but that we will 
ignore. The subtlety is related to the fact that the parton model (including the distribution 
function model of the hadronic structure of the photon) implicitly assumes that the con- 
stituents are massless (or, that the mass of the partons is small relative to all energy scales 
in the problem, and so can be safely ignored). As we are producing a massive 6-quark, the 
minimum s is of order ml, and it is not necessarily possible to ignore the initial 6-quark mass 
compared to all mass scales, namely Vs. In practice it turns out that our results are only 
sensitive to a finite rrib for two cases; at the lower energy colliders where the statistics are 
insufficient to make a meaningful measurement and for the low x region of da/dx distribu- 
tions (where, following the standard parton model notation, s = xys and pb = xyps, y is the 
fraction of the electron's momentum carried by the photon, and x is (nominally) the fraction 
of the photon's momentum carried by the 6-quark). In the latter case the uncertainties are 
only in the lowest x region where the statistical errors will still overwhelm the uncertainties 
in our definition of the scale. 

In the figures that follow, only the estimated statistical uncertainties are included in 
the error bars, so that the error bars shown will slightly underestimate the actual errors. 

A physically measurable quantity is the ratio of the initial 6-quark energy or momentum 

9 



to the beam energy. This will generally be closely related to r = xy, and the conversion 
is easily calculated. We show, in Fig. 6a, the number of events vs. t for various photon 
distribution functions in the 67 case. The largest differences in the models are in the lowest 
r bins where the event numbers are largest. Considering only the lowest r bin and using the 
Duke-Owens results to estimate the error we obtain a statistical error of roughly 2.5%. Thus, 
using only the lowest r bin for the measurement, and the higher r bins as a normalization, 
it will be quite easy to distinguish between the various distribution functions. Fig. 6b shows 
the same distribution in the e+e^ case. Here, the conclusion is not so clear, due to relatively 
larger error bars, but it should be possible to distinguish between LAC, DO/DG and GRV, 
though it may not be possible to distinguish DO from DC 

Finally we show, in Fig. 6a (67) and in Fig. 6b (e+e^) the distribution of event numbers in 
X. This distribution will be useful if one can tag the electron which provides the Weizsacker- 
Williams photon (allowing a complete reconstruction of the 67 initial state, including y), 
or in the 67 case if one can deconvolute the r distribution knowing the parton level cross 
section and the laser backscattered photon spectrum. As with the r distribution, the large 
X bins can be used to fix the normalization of the distribution while the lowest x bin can be 
used to provide information on the 6-quark distribution in the photon. Similar conclusions 
can be reached here as well: in the 67 case, it should be possible to distinguish between the 
various distribution functions, while in the e^e~ case it should be possible to distinguish 
between LAC, DO/DG and GRV (though possibly not between DO and DC). 

IV. CONCLUSIONS 

In this paper we proposed using high p^ 6-quarks as a means of determining the heavy 
quark content of the photon. We have shown that the process e + 7— >e + 6 + X provides 
a clean method of extracting the quark content of the photon and can distinguish between 
the different models in the literature. The subprocess we examined was e + b ^ e + b where 
the beam electron undergoes a hard scatter from a resolved b in the photon, in essence 
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Rutherford scattering off a photon target. The signal for the process is quite distinct; all 
backgrounds can be eliminated by insisting that only the scattered 6-quark balanced against 
the beam electron be observed and imposing a minimum p^ cut on the b to eliminate the 
remaining direct and once resolved bb production. 

It is clear that this reaction will allow a good determination of the quark content of 
the photon. In the low x-regions of the distributions, where the models differ the most, 
typical uncertainties are roughly 2.5% for the 67 colliders. This is more than adequate to 
distinguish between different models. For the e^e^ cases the statistics are poorer and do 
not offer quite as good a measurement, though valuable information can be extracted from 
these data. One could also use this approach to study the charmed quark content of the 
photon |2^ . The c-quark content of the photon is roughly four times as large as the 6-quark 
content due to the different quark electric charges. However, since charmed mesons are more 
difficult to reconstruct, the 6-quark gives a cleaner signal. The process we have proposed 
is complementary to others that have been studied in the literature and can contribute to 
more precise determinations of quark and gluon distributions inside the photon. Along with 
other processes considered elsewhere it is clear that an 67 collider is an ideal facility for 
studying the hadronic content of the photon. 
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FIGURES 

FIG. 1. Feynman diagram for the process 67 -^ eb + X. 

FIG. 2. Comparison of direct e + 7^e + 6 + 6 (dashed curve), e + [g]^ ^ e + b + b (dotdashed 
curve) and the sum (dotted curve) to the signal process e + 7 -^ e + b + X (sohd curve). Ah 
curves correspond to a 500 GeV e"'"e~ cohider operating in 67 mode. For the signal, e\b]^ -^ eb, 
(solid curve) we require that the final state electron and 6-quark be seen, 6^ more than 10" from 
the beampipe and | cos 0^1 < 0.85; for the direct processes, there is an additional requirement, 
|cos%| >0.85. 

FIG. 3. Diagrams for background processes, (a) Direct contributions; (b) Singly resolved con- 
tributions; (c) Doubly resolved contributions. 

FIG. 4. p^ distribution for the signal (solid line) and the most dangerous backgrounds: 
7 + [9']7 -^ b + b (rightmost dotted line), 7 + 7^6 + 6 (rightmost dashed line), [s']7 + [s']7 -^ b + b 
(leftmost dotted line), [9]7 + [^]7 -^ b+b (leftmost dashed line), e^ + e" -^ 7 + 6+6 (dotdashed line). 
These are all make use of the Duke and Owens photon distribution functions (where appropriate) , 
and are for a 500 GeV e+e^ collider operating in 67 mode. 

FIG. 5. Pj, distribution for the signal with LAC (solid line), GRV (dashed line), DO (dotted 
line) and DC (dotdashed line) distributions, using Q^ = s at (a) a 500 GeV e~^e~ collider and (b) 
a 500 GeV e+e^ collider configured as an 67 collider. 

FIG. 6. Event number vs. r for a 500 GeV e^e~ collider assuming L=50 fb^^ operating in (a) 
67 mode and (b) e+e^ mode, using the four different sets of photon distribution functions. The 
solid curve is for LAC set 1, the dashed curve for GRV, the dotted curve for DO and the dotdashed 
curve for DC. 
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FIG. 7. Event number vs. x for a 500 GeV e^e~ collider assuming L=50 fb^^ operating in (a) 
67 mode and (b) e^e~ mode, using the four different sets of photon distribution functions. The 
solid curve is for LAC set 1, the dashed curve for GRV, the dotted curve for DO and the dotdashed 
curve for DG. 
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